Competing orderings in an extended Falicov-Kimball model 
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We present a Hartree-Fock study of the Falicov-Kimball model extended by both on-site and non- 
local hybridization. We examine the interplay between excitonic effects and the charge-density wave 
(CDW) instability known to exist at zero hybridization, ft is found that the CDW state remains 
stable in the presence of finite hybridization; for on-site hybridization the Coulomb interaction 
nevertheless strongly enhances the excitonic average above its value in the non-interacting system. 
In contrast, for non-local hybridization, we observe no such enhancement of the excitonic average 
or a spontaneous on-site hybridization potential. Instead, we find only a significant suppression of 
the excitonic correlations in the CDW state. A phenomenological Ginzburg-Landau analysis is also 
provided to understand the interplay. 

PACS numbers: 71.30.+h, 71.28.+d 



The Falicov-Kimball model (FKM) describes a tight- 
binding system of itinerant d-electrons interacting via 
on-site Coulomb repulsion U with localized /-electrons 
of energy e/. The FKM was originally introduced as 
a minimal model of valence transitions in systems such 
as SmB 6 and Ce: by varying the inter-orbtial Coulomb 
repulsion U or the /-level e/, both discontinuous and 
continuous changes in the distribution of the electrons 
across the localized and itinerant states were foundi It 
was soon realised, however, that some overlap between 
the d- and /-wavefunctions was an essential feature of 
most systems displaying valence instabilities^ This "mix- 
ing" of the electron wavefunctions may be explicitly in- 
troduced by the inclusion of a hybridization potential V. 
A variety of methods, including Hartree-Fock^ real-space 
renormalization groups and alloy-analog approximation, 5 
revealed that the hybridization removed the previously 
observed discontinuous valence transitions. Work on the 
FKM ceased in the mid-1980s as it became apparent that 
the Periodic Anderson Model offered a more realistic de- 
scription of valence transition physics^ 

As interest in the FKM as a model of valence transi- 
tions waned, it was adopted as a model of a simple bi- 
nary alloyi In the limit of vanishing hybridization the /- 
electron occupation at each site is a good quantum num- 
ber: fixing the /- and d- populations, the ground state is 
identified as the configuration adopted by the /-electrons 
that minimizes the energy of the conduction electrons. In 
particular, for a bipartite lattice at half-filling and equal 
concentration of d- and /-electrons, the /-electrons oc- 
cupy the sites of one sublattice only, the so-called che- 
qucrboard phase. For dimension d > 2, this chequer- 
board charge-density wave (CDW) state obtains for tem- 
peratures below a critical temperature Tqdw] above this 
temperature a disordered phase is realized. For d = 1 the 
critical temperature is zero. We note that the FKM as 
a binary alloy has been extensively studied in the case 
of infinite dimension d — > oo: the dynamical mean-field 
theory (DMFT) gives an exact solution in this limit £ 



The FKM with hybridization has lately attracted re- 
newed attention due to the investigation of optical prop- 
erties in this model by Portengen et almi Following closely 
Leder's Hartree-Fock (HF) work^ they found that the 
Coulomb repulsion induced an effective on-site hybridiza- 
tion; this effect was sufficiently strong that it persists in 
the limit of negligible hybridization. In fact, their cal- 
culations were performed exclusively in this limit: their 
solution with non-zero polarization or excitonic average 
(d\f) is indistinguishable from the well-known excitonic 
insulator (EI) stated The "spontaneous" excitonic av- 
erage was interpreted as evidence of electronic ferroelec- 
tricity. Their HF solution, however, assumed a homoge- 
neous ground state for the system; the possibility of a 
CDW ground state was not considered. 

The problem of reconciling the results of Portengen 
et al. with the known CDW instability has only been 
partially addressed. Since the DMFT equations are no 
longer exactly solvable for non-zero hybridization poten- 
tial, CzycholUi performed a HF analysis for the d — > oo 
model. It was found that for V = there was no spon- 
taneous excitonic average and that the CDW phase was 
stable against sufficiently small on-site hybridization. For 
a given U there was a critical hybridization V C (U) un- 
der which the CDW phase prevails. Czycholl neverthe- 
less concluded that the the inter-orbtial U could strongly 
renormalize the hybridization, and so could be important 
in the description of the optical properties of strongly 
correlated electron systems. Also working in the limit 
of large spatial dimensions, Zlatic et al. noted that 
for V = the hybridization susceptibility diverges as 
T —y 0, although they concluded that a generalization 
of the FKM would be required for (d* f) ^ at finite 
temperatures^ 

Comparatively little work has been done on this prob- 
lem in finite dimensions. Farkasovsky has used exact- 
diagonalization and the density matrix renormalization 
group methods on small one-dimensional systems to rule 
out the possibility of a spontaneous excitonic average at 
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zero temperaturei* By the same methods, Farkasovsky 
has also analyzed the effect of locaiii and non-locali^ 
hybridization; these works are more concerned with the 
effect of the hybridization on valence transitions, ignor- 
ing the possibility of an excitonic renormalization of the 
hybridization potentials. Batista and co-workers have 
claimed that a non-local hybridization stabilizes ferro- 
electricity in a FKM extended by /-hoppingfi^ Sarasua 
and Contincntino have investigated a similar system^ 

The FKM extended by hybridization cannot be solved 
exactly and so it is necessary to use approximate methods 
to understand the properties of the model. In this paper, 
we present a HF study of the effect of the hybridization 
upon the CDW state on a two-dimensional square lat- 
tice. The HF approximation is reliable for small temper- 
atures. It tends, however, to overestimate the stability 
of ordered phases: in particular, the HF result for crit- 
ical temperature Tqdw is very likely to be larger than 
the exact value. Nevertheless, we can reasonably expect 
that the HF approximation will give at least a qualita- 
tively correct account of the relative stability of ordered 
phases, even in 2D. The HF is therefore an appropriate 
tool to study the competition between the EI and CDW 
phases in the FKM. We consider only e/ = and half- 
filling (the particle-hole symmetry point) as the CDW 
state here adopts the simple chequcrboard form; for these 
parameters also the excitonic average takes its maximum 
as shown in the analysis of Portengen et al& 

The FKM Hamiltonian for spinless Fermions is written 

H = -tJ2 4d 3 +e f J2 n f j + £{^4^ + H.c.} 



(i,3) 



external magnetic field. In this case, we cannot speak 
of an EI phase, but rather an excitonic enhancement of 
the hybridization. This will be apparent if A exceeds its 
value in the U = system. The modulation factor Aq 
is included in Eq. for completeness. In the usual HF 
treatment^ a homogeneous solution is assumed and so 
5d = Sf = Aq = for all values of the Coulomb interac- 
tion. 

We thus obtain for the HF Hamiltonian 



H HF = -tJ24 d i ■ +£/^(n / + <5/cos(Q 

(id) 3 



+Uj2( nd + S d cos(Q-r J ))r 



E{(^' -U[A + A q cos(Q • rj)] SiMfj 



i j 
-H.c 



(5) 



An important feature of this Hamiltonian is the mean- 
field renormalization of the d-f hybridization poten- 
tial by the inter-orbital Coulomb interaction, Vij — > 
Vij — U[A + Aq cos(Q • rj)] Sjj. The effective on-site hy- 
bridization potential introduced by the decoupling of the 
interaction is responsible for the spontaneous polariza- 
tion in Portengen et aVs work. 7i HF is diagonalized by 
the canonical transform 



(1) 7 r = <d k + <d k+Q + sr/k + cr/k+q , (6) 



Some overlap between the d- and /-electron wave- 
functions is assumed, hence the hybridization term 



The concentration of electrons is fixed at 1 



7f ) wncrc N is the number of sites. We 

measure all energies in terms of the d-electron hopping 
integral t. 

In our HF decoupling of the Coulomb interaction, wc 
include the possibility of the CDW state by allowing for 
a periodic modulation of the order parameters: 



/ C0S(Q • Tj) 



(nj) = rJ 
{nf} = n d + <5 d cos(Q • Yj 
(fUj) = A + A q cos(Qt 



(2) 
(3) 
(4) 



The nesting vector Q = (-, — ) where a is the lattice 
constant. The order parameter of the CDW state is Sd 
and 5f for the d- and /-electrons respectively. Note that 
we require sgn(£/) = — sgn(Sd)- A is the excitonic aver- 
age; in the absence of an on-site hybridization potential 
V, A ^ indicates the EI phase. When V ^ 0, the 
El-normal phase transition is lifted from criticality, in 
analogy to the fcrromagnct-paramagnct transition in an 



where m = 1,2,3,4. The coefficients in Eq. (jBJ) are ob- 
tained by solving the associated Bogoliubov-de Gcnncs 
(BdG) eigenequations: 



(7) 



where 



and 



/ e k + Un f US f Vk-C/A ~UA Q \ 
US f e k+Q + Unf -UA Q V k+Q -UA 
UA* Q Un d US d 

US d Un d J 

(8) 



-ua* q y k * +Q - II A* 



m 1} rn £ m /-m\Transpose 



<x f m / m m cm /-m\ 

*k = ( w k , " k , £k , Ck ) 



(9) 



Here £k = — 2t(cos(k x a) + cos(k y a)) is the d-electron en- 
ergy dispersion. The self-consistency equations for the 
HF parameters may be written in terms of the BdG eigen- 
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vectors: 
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N 
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k E ' E + /(^) ■ (io) 

k m 

^E'{(4 + q4) + (44 + q}} 
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k m 
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k 

^ E ' E <cr + srcn /(^n • (13) 

k m 
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k 

= ^E'E^r<+cr«n/(^D- (14) 

k m 

A « = ^E'{(/k + Q rf k) + (/k^ +Q )} 

k 

= ^E'E^^+ero/^D • (15) 



A = 



k m 



The prime denotes summation over half the Brillouin 
zone; f(E) = 1/{1 + cxp[/3(£' — /i)]} is the Fermi dis- 
tribution function. The chemical potential fi is deter- 
mined by the condition 1 = iV" 1 £ k ' £ m We 
use an exact diagonalization method to solve the BdG 
equation 10 self-consistently. We start with an initial set 
of order parameters. By solving Eq. (J7J), the new order 
parameters are computed via Eqs. I|1U|) to (|15|) and are 
substituted back into Eq. Q). The iteration is repeated 
until a desired accuracy is achieved. 

We first consider the case of an on-site hybridization, 
Vij = VSij . In agreement with previous workii*i2ii^ we 
find that for vanishing hybridization the CDW phase is 
always stable against the EI phase and there is no spon- 
taneous excitonic average. The CDW order displayed by 
the V = ground state will persist in the presence of 
sufficiently small hybridization potentials, although the 
transition temperature Tqdw will t> e considerably sup- 
pressed, see Fig. ^ We find, however, that for finite 
hybridization the Coulomb interaction will strongly en- 
hance the magnitude of A. We plot the variation of |A| 
with temperature in Fig. [2 Note that since A = — |A| 
there is a large renormalization of the hybridization po- 
tential due to the mean-field decoupling of the Coulomb 
interaction Eq. (j3J). Comparing the homogeneous so- 
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FIG. 1: Variation of the CDW order parameter Sd with tem- 
perature in the absence (solid line) and presence (dotted line) 
of an on-site hybridization potential. 
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FIG. 2: Comparison of the excitonic average |A| in the ab- 
sence (solid line) and presence (dotted line) of the CDW in- 
stability for U = l.Ot with the value in the non-interacting 
system (dashed line) for V = 0.2t. 



lution without the CDW ordering (the solid line) with 
the solution with a coexisting CDW ordering (the dot- 
ted line), we find a significant suppression of |A| at the 
onset of the CDW order at T ~ O.lt. Even within the 
CDW phase, however, the excitonic enhancement of the 
on-site hybridization is still apparent as |A| exceeds its 
value within the non-interacting system (the dashed line 
in Fig- Ell ■ We do not find any evidence of non-zero Aq. 

This competition can be understood from a phe- 
nomenological Ginzburg-Landau (GL) theory. The GL 
free energy density, in terms of both the CDW (Sd) and 
EI (A) order parameters, can be constructed from a sym- 
metry analysis: 



a E i\A 
+P3\A\ 2 \5 d \ 2 



acDw\8d\ 2 



/3 4 (A* 2 <5 2 



■ falAf + f3 2 \6 d \ 2 
A 2 5f) , (16) 
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where we assume clei = o-'ei^ 
a CDw( T ~ t cdw)- We assume , 



T EI ) and a C DW 
(i = 1,2,3,4) are all 
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positive. In the region where T EI > T CDW ^ the second 
phase transition temperature for the CDW ordering is 
renormalized by the pre-existing EI order parameter: 



Tcdw 



1 CDW 



(/3 3 -2/3 4 )(Tg / 



1 CDW) 



Wia' CDW l<x' EI - (As - 2/3 4 ) 



(17) 



It means that when the EI order parameter pre-exists, 
the second phase transition temperature for the appear- 
ance of the CDW order parameter can be strongly sup- 
pressed by the dominant EI order parameter. This ex- 
plains why the transition temperature Tcdw decreases 
with increased hybridization potential V, as shown in 
Fig. ^ Below the second phase transition temperature 
Tcdw , a little algebra yields 



A 
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-2f3 2 UEI + ucDwifc - 204 
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4/3i/? 2 - (Ps 



2/3 4 ) 2 



1/2 



1/2 



(18) 
(19) 



Under the condition that the temperature derivative 
ct'cDW IS l ar g er than a' EI , which is indeed confirmed by 
our numerical results near Tcdw (see Figs. 11121) . a CDW 
changes more rapidly than uei when the temperature is 
lowered. Consequently, the CDW order 5d increases while 
the EI order A decreases with the lowered temperature. 

Despite the popularity of the on-site hybridization po- 
tential, this is actually forbidden in real d-f systems by 
parity considerations^ We are instead required to con- 
sider a non-local hybridization with inversion symmetry: 
the simplest such potential is 



'<lf 



l y3y ' 



l yJy 



ft. 



-S. 



(20) 



where any site on the lattice is given by = i x ax. + i y ay. 
This is a particularly interesting case as the Coulomb- 
induced hybridization has a different (s-wave) symmetry. 
In the non-interacting system, the (on-site) excitonic av- 
erage A vanishes; the non-local hybridization potential 
instead gives rise to an anisotropic excitonic average 



sin(k x a) + sin(k y a)) (/^k) \ ■ (21) 



The study of this quantity allows us to assess the effect 
of the inter-orbital Coulomb repulsion upon the d-f hy- 
bridization. 

As with the on-site hybridization, we find that for given 
U the CDW phase is suppressed by the presence of the 
non-local hybridization (see Fig. |3J). We do not, how- 
ever, find any evidence for a Coulomb-induced on-site 
hybridization when the CDW instability is allowed: for 




FIG. 3: Variation of the CDW order parameter Sd in the 
absence (solid line) and presence (dotted line) of a non-local 
hybridization potential. 
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FIG. 4: Comparison of the anisotropic exciton average H in 
the non-interacting and interacting systems. We have tdj = 
0.2t. 



all non-zero tdf we have A = Aq = 0. Czycholl con- 
sidered the appearance of an on-site average A to be 
likely due to the substantial excitonic enhancement of the 
on-site hybridization potential by the Coulomb interac- 
tioniii Our results clearly demonstrate that this El-like 
scenario, and the consequent formation of an electronic 
ferroelectric state, is severely compromised by the pres- 
ence of non-local hybridization. 

In Fig. 0] we plot S as function of temperature in both 
the interacting and non-interacting systems. The onset of 
CDW order for the given Coulomb values occurs at the 
point of intersection of the broken lines with the non- 
interacting (solid) line. Remarkably, for the standard 
homogeneous solution there is no effect on S due to the 
Coulomb interaction: the variation of S with temper- 
ature exactly follows the curve for the non-interacting 
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system. Within the CDW phase, however, S is sup- 
pressed below its value in the non-interacting system. 
We offer the following explanation for this anomaly: 
the hybridization potential Eq. I|2U|) connects the A- 
sublattice <i-orbitals with B-sublattice /-orbitals and vice 
versa. Assume that in the CDW state the A-sublattice d- 
orbitals have n d > 0.5 and so the B-sublattice /-orbitals 
have n' > 0.5; clearly A-B sublattice d-f hopping will 
be suppressed, hence also the reduction in S. 

In conclusion, we have examined the competition be- 
tween excitonic and CDW instabilities in the FKM ex- 
tended by both on-site and non-local hybridization. In 
both cases, we find that the CDW phase remains sta- 
ble at low temperatures even in the presence of a fi- 
nite hybridization. For the local hybridization we find 
that the Coulomb interaction nevertheless strongly renor- 
malizes the hybridization potential in agreement with 
previous work^i The situation is qualitatively different 
for the more realistic non-local hybridization: there is 



no enhancement of the non-local hybridization and the 
Coulomb interaction does not induce a spontaneous on- 
site hybridization. Within the CDW phase, the non-local 
hybridization is suppressed in line with the increasing lo- 
calization of the d- and /-electrons. The failure of the 
Coulomb interaction to induce an effective on-site hy- 
bridization except when such a term is already present 
casts significant doubt over the usefulness of Eq. Q as 
a minimal model for electronic ferroelectricity. Inter- 
orbital Coulomb repulsion may nevertheless still be im- 
portant for understanding optical properties of strongly- 
correlated electron systems: for example, a recent exten- 
sion of the FKM by /-electron hopping offers a plausible 
scenario where the formation of an exciton BEC gives a 
spontaneous excitonic averaged 
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